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Abstract 

A decay of the gravitino dark matter is an attractive candidate to explain the current excesses of 
the PAMELA/ATIC cosmic-ray data. However, R-parity violations are required to be very tiny 
in low-energy scale. We suggest a R-parity violation in the right-handed neutrino sector. The 
violation is suppressed by a see-saw mechanism. Although a reheating temperature is constrained 
from above, the thermal leptogenesis is found to work successfully with a help of the R-parity 
violating right-handed neutrino. 
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I. INTRODUCTION 



Observations of comic rays have been greatly developed. After a series of cosmic-ray 
measurements, PAMELA recently published the first result of the positro^ fraction, which 
show an excess compared to the background at energies above lOGeV [1]. Interestingly, 
the ATIC collaboration also reported an excess of the electron plus position flux in a range 
of O(100)GeV [2]. Although the experimental data and the background estimations still 
involve ambiguities, the excesses may be a sign of the dark matter (DM). In this letter, we 
focus on decaying DM scenarios. The DM is not always required to be stable as long as its 
life time is long enough for the DM to survive until today. Then, its decay products are a 
possible source of the energetic cosmic rays. 

It is non-trivial to realize such a long-lived particle. Actually, decay operators must 
be incredibly suppressed for the models to be viable. A well-motivated candidate of the 
DM is known to be the gravitino with a broken R-parity. This is because its decay is 
doubly suppressed by the Planck scale and R-parity violations [3]. Even with the Planck 
suppression, the R-parity violation is limited to be very tiny and looks unnatural unless some 
mechanisms are taken into account. As a solution of the fine-tuning, we propose a R-parity 
violation introduced in the right-handed neutrino sector. Then, the violation is suppressed 
in the low-energy scale by a see-saw mechanism. Actually, since the violation appears in 
low-energy phenomena though the right-handed neutrino, its effect becomes suppressed by 
the right-handed neutrino mass scale. We will show that the PAMELA/ ATIC excesses are 
explained for M^v ~ lO^GeV without a parameter tuning. 

This type of the R-parity violation does not only explain the PAMELA/ATIC anomalies, 
but also opens a window for the successful thermal leptogenesis j4| especially in a large 
gravitino mass region. Indeed, the mechanism is an attractive and elegant solution to the 
mystery of the baryon asymmetry of the universe, though it usually requires a relatively high 
reheating temperature, lO^GeV (see e.g. [5|])- In such a high reheating temperature, 
the severest cosmological constraints are provided by the big-bang nucleosynthesis (BBN) 
and the overclosure of the universe. The former constraint is due to a prolonged decay of 
the next-to-lightest superparticle (NLSP) into the gravitino lightest superparticle (LSP). In 
contrast, the NLSP can decay into the Standard Model (SM) particles much more efficiently 
through the R-parity violating operator. Thus, the BBN bound becomes ameliorated j^. 
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The overclosure constraint then puts a constraint on the reheating temperature. In par- 
ticular, the thermal gravitino production is enhanced in a larger soft mass region. According 
to PAMELA/ATIC, the cosmic-ray excesses imply a heavy gravitino, namely large super- 
particle masses. We will see that the leptogenesis temperature is unlikely to be satisfied 
for m3/2 ^ ITeV unless the gluino mass parameter is relatively small at the GUT scale. 
Nevertheless, since the R-parity violation in this letter contributes to a CP-asymmetric 
right-handed neutrino decay the thermal leptogenesis will be shown to work to explain 
the observed baryon asymmetry easily. Thus, we will study the R-parity violating scenario 
in the right-handed neutrino sector as an attractive framework of the decaying gravitino 
DM. 



II. R-PARITY VIOLATION IN RIGHT-HANDED NEUTRINO SECTOR 

Let us consider the gravitino DM scenario with a decay due to R-parity violations. A 
spectrum of the cosmic rays from its decay products are determined by the type of the 
R-parity violating operator, the gravitino mass and its life time. In a case of the bilinear 
R-parity violation, the dominant decay channels are V'3/2 — ^ Wi, Zu, 7Z/ and hu 
The decay is parametrized by the sneutrino vacuum expectation values (VEV). With the 
operator, W = fi'LH, it becomes 

{-) - (1) 

where L and H are the left-handed lepton and up-type Higgs, respectively. The parameters 
are the higgsino mass fi, the Higgs VEV v ~ 174GeV, the Higgs VEV ratio tan P = {H) / {H) , 
and the sneutrino soft mass Kjoft = rnli)*!). Then, the gravitino life time is estimated as 

for m3/2 ^ lOOGeV, where r] is defined as = \{i')\/v. According to the current 
PAMELA/ATIC data, the mass and life time of the gravitino is favored to be 

m3/2 > 0(100 - 1000)GeV, ~ lO^Sec. (3) 

which leads to the R-parity violation parameter as 
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This 7] satisfies with the other cosmological constraints [o^. The gravitino survives until 
today as 7] <^ iQ-i^-^) for 1113/2 = O(100)GeV, and the parameter is small enough to avoid 
the wash-out of the lepton/baryon asymmetry since rj < 10^"^ is satisfied. Furthermore, the 
BBN constraint is solved for rj^ iQ-i'^^-^'^) by virtue of the R-parity violation, and the R- 
parity violating contribution with rj < 10^^ is negligible for the neutrino mass. We assume 
that the neutrino mass is obtained by the see-saw mechanism in this letter. 

However, from the naturalness point of view, the violation is too tiny and must be fine- 
tuned. We propose a simple framework to realize this tiny R-parity violation. We introduce 
the following R-parity violating operator in the right-handed neutrino sector, 

W = XNHH. (5) 

Accompanied by the right-handed neutrino mass and the Yukawa terms, W = MiqNN/2 + 
Y]\fNLH, after decoupling the heavy right-handed neutrino, we obtain the effective R-parity 
violating operators in low-energy scale: 

w = _^^MHnmi (6) 

It is noticed that the R-parity violating effects appear with l/M^v because they contribute 
through the right-handed neutrino. Thus, the R-parity violating effects are naturally sup- 
pressed by the see-saw mechanism. 

Taking the Higgs VEVs, ([6]) behaves as a bilinear R-parity violating operator. Actually, 
combined with the /z term and the sneutrino soft mass, the sneutrino acquires VEV as 

. XYNfiv^sin^ (3 , . 

(u) ~ K . (7) 



Then, i] is estimated to be 

r? ^ 0.6 X 10"^° ■ A 



where m,y is a typical scale of the light neutrino masses defined as Ya? = V M]ym,j/ v sin {3. 
Here, we omitted flavor indices. As a result, we find that the gravitino life time T3/2 ~ lO^^sec 
required from PAMELA/ ATIC is realized for A ~ 1 and Mn ~ lO^GeV. We want to 
emphasize that there is no fine-tuning in the R-parity violating parameters. 

The gravitino decay is almost the same as that by the bilinear R-parity violation. In 
fact, it is dominantly induced by the sneutrino VEV ([7]), while the other decaying channels 



from ([6]) are subdominant. Although W = —X'^{HH)'^/2Mn is also derived from ([5]), the 
resultant decay operators are negligible for the cosmic ray spectra. 

According to the current cosmic-ray data, the positron spectrum has a steep rise in high- 
energy regiom This behavior prefers the gravitino DM decaying direct into the positron or 
anti-muon 
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12l |. The decay product is determined by flavour structure of the sneutrino 
VEV. Denoting the flavour indexes explicitly in ([7]), (uj) is proportional to \i{Y]^)ij/ Mj^i 
with respect to the flavor structure, where i is the index of the heavy neutrino and j for the 
light one. Thus, it is not surprising to expect the positron and anti-muon to be produced 
in the decay. Consequently, the R-parity violating operator ([5]) can naturally explain the 
current PAMELA/ATIC excesses with mg/a = 0(100 - 1000) GeV. 

Let us comment on radiative corrections with the R-parity violation They induce 
a linear term of the right-handed neutrino in the Kahler potential. This then leads to the 
right-handed sneutrino VEV in the framework of the supergravity. Although this can give 
a similar contribution to the effective R-parity violating operators ([6]), the results in this 
letter do not change qualitatively. Thus, we will neglect it hereafter. 

We may also have the R-parity violating operators in the right-handed neutrino sector 
other than ([5]) such as the superpotential terms, and A^^. However, their couplings are 
required to be suppressed since they can induce too large R-parity violations through the 
right-handed sneutrino VEV within the supergravity. Thus, we consider that the operators 
are forbidden by (discrete) symmetries. 

In the above, we focused on the SUSY invariant operators for the R-parity violations. 
Apart from them, we may have the violations in the SUSY breaking term. It can be noticed 
that the analysis is similar to the above, and the result is almost the same. 



III. UPPER BOUND ON REHEATING TEMPERATURE 

In the decaying gravitino DM scenario, the cosmological constraint comes from the over- 
closure of the universe. Actually, the gravitino is thermally produced in the hot plasma. At 
the leading order, the relic abundance is evaluated as 

ny,h ^ ^ c..,. (Tn) In ^ (^1 + j (yMg^ j ) ' 0) 
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where the definition of the parameters are found in [13|. In the numerical analysis, we 
include the electroweak contributions, which can be sizable especially when the gluino mass 
parameter is rather small compared to the Bino and Wino ones at the GUT (namely Tr) 
scale. It should be mentioned that the gravitino production rate includes an 0(1) uncertainty 
from unknown higher order contributions and nonperturbative effects pj,]. In addition, 
resummation of thermal masses potentially increases the rate by about a factor of two 

n 

[Iq . Since the massive gravitino behaves as a cold dark matter, its abundance is constrained 
from the measurements. According to the WMAP 5-year data, the abundance is required 



to satisfy [l7l |. 



^3/2/1^ < ^Duh^ ^ 0.1223, (10) 

at the 2a level. Thus, provided the gaugino and gravitino masses, we obtain an upper bound 
on the reheating temperature. 

It is interesting to study the overclosure bound in terms of the superparticle masses at 
the weak scale, which are expected to be measured by LHC From (Q we notice that the 
gravitino production is evaluated with the parameters at the reheating temperature scale. 
Then, they are correlated with those at the weak scale by solving the renormalization group 
equations. Let us mention that the 1-loop result includes large uncertainties especially in the 
colored sector: 2-loop contributions can give an O(10)% correction, and the renormalization 
scale at the 1-loop level potentially contains additional O(10)% uncertinty, which is reduced 
by taking the higher order corrections into account. Noting that (j9]) includes the gluino 
mass squared, the thermal gravitino abundance can change drastically. Actually, 1)3/2/1^ is 
found to be almost as twice as the previous 1-loop analyses for Mi{T^) ^ ^3/2 with the 
gluino mass fixed at the weak scale [l^ (see also the discussion in [18]). In the following 
numerical analysis, we evolve the renormalization group running at the 2-loop accuracy and 
include the 1-loop threshold corrections by means of SOFTSUSY 2.0.18 jl^. 

In Fig. [H we plot the maximal reheating temperature allowed by the overclosure bound 
with varying the gluino mass M3 compared to the gaugino mass M1/2 at the GUT scale. 
Here, the Bino and Wino masses are set to be equal Mi = M2 = M1/2 and realize the 
NLSP mass mNLSP on each solid lines. Since the universal soft scalar mass mg except for the 



Apart from the thermal production, we also have nonthermal contributions to the gravitino abundance, 



particularly from inflaton decay 16|. Since they are model dependent, we neglect them for simplicity. 
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Higgs ones is chosen to be mo = M1/2, the Bino-hke neutrahno is the NLSP, while the Higgs 
mass squareds are assumed to be = mg and = —ttlq for the electroweak symmetry 
breaking to give arise especially in a large mg region. The scalar trilinear couplings, Oq, 
and tan/5 are = and tan/5 = 30, though they are irrelevant for the maximal reheating 
temperature. 

The gravitino mass is relevant for the upper bound on the reheating temperature. We 
assumed that it is equal to the NLSP mass. This turns out to maximize Tr for M3/M1/2 > 
0.5. Although the possible maximal temperature is obtained by a smaller gravitino mass for 
a lighter gluino region, the difference of Tr is less than 10% even for M3/M1/2 = 0.2. It is 
commented that the gluino can become the NLSP for M3/M1/2 < 0.2. 

We find that the maximal reheating temperature can be as large as O(10^)GeV in a 
large parameter region, while Tr > lO^GeV is realized for a smaller NLSP mass, e.g. 
"^NLSP ~ SOOGeV for the universal gaugino mass M3 = M1/2. This is because the ther- 
mal gravitino abundance ([9]) increases as the soft mass scale is larger. From Fig. [H we 

I for a hier- 



see that Tr > lO^GeV leads to mNLSP < l-STeV for M3/M1/2 > 0.2 (see 
archical gaugino mass case). According to PAMELA, the gravitino mass is indicated as 
"^3/2 ~ O(100)GeV, i.e. rriNLSP ~ O(100)GeV. Then, it is easy to obtain the leptogenesis 
temperature in a wide class of models (see [3], in which discussions of superparticle mass 
spectrum are also explored). 

In contrast, the ATIC result implies a larger gravitino mass, m3/2 ^ ITeV. Then, the 
gluino mass is required to be suppressed. For instance, when the gravitino mass is ITeV, 
Tr > lO^GeV needs M3 < 0.7Mi/2 at the GUT scale, and M3 < 0.3Mi/2 for 7713/2 = 1.5TeV. 
Thus, the SUSY breaking models become specified to realize the leptogenesis temperature. 
A lighter gluino mass parameter is favored at the GUT scale, and thus a mass spectrum of 
the superparticles tends to degenerate at the weak scale especially when the gravitino mass 
is heavier ^. On the other hand, the universal gaugino mass models look inconsistent with 
the thermal leptogenesis in the light of the ATIC result, though the usual gauge- and many 
gravity-mediation models have this feature. Thus, we will revisit the thermal leptogenesis 
in the next section in the presence of the R-parity violation. 



^ This mass spectrum looks like that of the mirage mediation 26[. However, it is usually obtained by 
enhancing the anomaly mediation, namely with a large gravitino mass, e.g. lOOTeV. Thus, the 

gravitino is not the LSP. 
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FIG. 1: The maximal reheating temperatures allowed by the overclosure bound are drawn for 
fixed NLSP masses. The gravitino mass is set to be equal to the NLSP mass. Note that these 
masses are the physical one, while M3 is given at the GUT scale, and the soft scalar masses are 
chosen to be mo = M1/2 with = Mi = M2 at the GUT scale. Then, the Bino-like neutralino 
is the NLSP for M3 ^ 0.2Mi/2- See the text for the other irrelevant parameters. 

Before proceeding to the next section, let us comment on the other parameter dependence. 
In the analysis, we take mo = M1/2 at the GUT scale. For tiny mo, the stau can be the 
NLSP. Then, the upper constraint in Fig. [1] becomes severer for fixed mNLSP? since the 
gaugino masses increase. Thus, the model parameters discussed above are limited more 
severely. 

IV. THERMAL LEPTOGENESIS REVISITED 

We now revisit the thermal leptogenesis in the presence of the R-parity violating operator 
In the above, we found that a gluino mass is bounded from above to realize the leptoge- 
nesis temperature Tr ^ lO^GeV in a large gravitino mass region. In particular, the universal 
gaugino mass models are inconsistent with the temperature if we take the ATIC result. 
On the other hand, we have introduced the R-parity violating operator in the right-handed 
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neutrino sector. We will see that a CP asymmetry of the decay of the right-handed neutrino 
is enhanced by the violation Q, and thus the leptogenesis mechanism works successfully. 
In the thermal leptogenesis, the baryon density relative to the photon density becomes 

(cf. Q) 

— ~ -1.04 X lO'^eiK, (11) 



which is compared to the observation 



id, 



— = (6.21 ±0.16) X 10"^°. (12) 

Here and hereafter, we assume the hierarchical right-handed neutrinos and ignore the flavour 
effects which does not affect the lower bound on the reheating temperature. The efficiency 
factor K represents effects of washout and scattering processes, which is obtained by solv- 
ing the Boltzmann equations. For the case of zero initial abundance of the right-handed 

nn 

neutrinos, its maximal value is known to be k ~ 0.2 [5|, |21[. 

The CP asymmetry ei in the decay of the lightest right-handed neutrino, A^^i —>■ LH, 
is determined by the structure of the neutrino sector. With k^O.2 the observed baryon 
asymmetry requires 

|ei| > 3 X 10"^. (13) 



In the R-parity preserved models, it is known that the |ei| has an upper bound as 22,12 



r(A^i ^ L + H)- T{Ni ~^L^ + H^ 



< 



87r(iJ)^ mi + ms 



(14) 



T{Ni L + H) + T{Ni ^ L^ + H^ 

due to the see-saw relation. Here mj (mi < m2 < m^) are the mass eigenvalues of the 
light neutrinos, and M^i is the mass of the lightest right-handed neutrino A'^i. Using the 
atmospheric neutrino mass squared difference Am^^j^ ~ (2.5 ± 0.2) x 10~^eV^, the lower 
bound on the right-handed neutrino mass becomes 

Mmi > 1.4 X 10^ GeV ■ sin^/? (15) 

with the 3(7 values of nsln^ and Am^^.^^. Consequently, the corresponding lower bound on 
the reheating temperature is obtained as Tr ^ 1 x 10^ GeV. 

The upper bound on |ei| f fM]) can be relaxed by the R-parity violating operator ([5]). The 
operator contributes to the CP asymmetric decay, A'^i LH, via the diagrams exchanging 
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H and H (and A^) in the loop. In addition, a decay, A^i — > HH, gives arise with the R-parity 
violating operator of Ni. Consequently, we obatin j?! 

1 1 



.(RPV) ^ 

8^(l^jvl^i^)ii + |Ai| 



X Im 



(16) 



where the loop function is /(x) = — y^ln(l + l/x) — 2A/x/(a; — 1). We notice that the result 
depends on Ai and Aj (i 7^ 1). In order to avoid a strong washout by a (inverse) decay, Ai 
is favored to satisfy |Aip ^(YAr^iv)ii' while Aj is allowed to be 0(1). Then, the R-parity 
violating processes can dominate the CP asymmetry of the right-handed neutrino decay. 
Actually, it is estimated as 



e^''^ ^ 2 X 10-^ ■ Im[ctA,, 



1 _ 1 _i 

2 



m' 



(17) 



VlOSCeVy V 10 y \10-^eV J VO-leVy sin^/?' 

for M]\ri/M]\ri ^ 1. Here, the parameters are defined as (Y]\rY^)ii = rfi'^yj MniMn-i/v"^ sin^ /3, 
mi = [(Fat F^) 11 + \ \i\^]v'^ / M^i and Ci = XiV / ^/Wmml. We can see that the result exceeds 
the upper bound (|T4l) for Aj ~ 1 and |Aip ^(yAry^)ii, i.e. ci ~ 1. 

In the last section, we found that Tr = O(10^)GeV is obtained even in a heavy NLSP 
mass region from Fig. [TJ On the other hand, we saw in this section that the CP asymmetry 
of the right-handed neutrino can be drastically enhanced by the R-parity violating operator 
([5]), and the condition (|T3l) is easily satisfied, e.g. for Mi = lO'^GeV. Since A''i is expected to 
be produced sufficiently in the thermal bath once Tr exceeds Mati, the thermal leptogenesis 
consequently works in a wide class of models with explaining the cosmic-ray anomalies from 
PAMELA/ATIC. To be explicit, let us give an example: Aj^i ~ 1 with the lightest right- 
handed neutrino mass Mni = O(10*)GeV and Mni/Mni = 0(10). 

One may be worried that such a large R-parity violation enhances washout effects. The 
(inverse) decay and AL = 1 processes indeed receive additional contributions from the 
R-parity violating operator. Considering that they are controlled by the washout mass 
parameter, mi = {YmYIi)iiv'^ / M^i, we notice that the extra effects just redefine mi as the 
new one given in the previous paragraph. On the other hand, the R-parity violation can 
enhance the AL = 2 washout. Actually, we obtain LH HH {H means H or H) by 
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exchanging a heavier right-handed neutrino in the diagram. Since the diagram is almost 
the same as the R-parity preserved one except for a couphng and corresponding field, the 
additional contribution is roughly estimated as 

(sub) \-^i\'^\O^N)il\'^ Mjf^ . . 



where 7^2=2 represents the AL = 2 washout term with a subtraction of the resonance 21 1. 
Thus, 7rpv becomes larger than ^^^^=2 f*^^ ^ lO^^GeV. Since 7al=2 negligibly small 
for this Mni, we can numerically check that the washout is dominated by the (inverse) decay 
and AL = 1 scattering even for Aj = 0(1) (see e.g. 2l| for the numerical estimations of 
the washout terms in the R-parity preserved case). Although 7^2=2 increase as M^i grows, 
7rpv decreases at the same time. Thus, the R-parity violating contribution does not change 
the result. 



V. CONCLUSIONS 

The gravitino DM scenario is an attractive candidate to explain the PAMELA/ATIC ex- 
cesses. However, the R-parity violation has to be tightly suppressed in the low-energy scale. 
In this letter, we proposed a framework of the R-parity violation to solve the fine-tuning. 
Introducing the violation in the right-handed neutrino sector, the effects are naturally sup- 
pressed in the low-energy scale by the see-saw mechanism. We found that the life time 
required by PAMELA/ATIC, T3/2 = O(10^^)sec, is obtained by the R-parity violating oper- 
ator W = XNHH with A ~ 1 for = O{10^)GeV. 

The R-parity violation not only provides a well source of the cosmic rays to explain the 
PAMELA/ATIC anomalies, but also opens a window for the thermal leptogenesis. In the 
presence of the R-parity violation, the reheating temperature is bounded from above by the 
overclosure bound due to the thermal gravitino production. We found that Tr ^ lO^GeV 
is accessible for 712^/2 < 600GeV, while ^ ITeV favored by ATIC restricts the gluino 
mass to be suppressed for Tr ^ lO^GeV. Nonetheless, the leptogenesis was shown to work, 
since the R-parity violating interactions in the right-handed neutrino sector enhances the CP 
asymmetric right-handed decay. Therefore, the gravitino LSP with the R-parity violation 
in the right-handed neutrino sector is a nice candidate of the decaying DM scenarios. 

It is worthwhile to clarify the parameter dependence of the phenomena. The cosmic-ray 
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spectra are sensitive to larger A2 and corresponding M7V2 but insensitive to Mi^i as long as 
Ai ^ A2, while the lightest right-handed neutrino determines the thermal leptogenesis, i.e. 
the baryon asymmetry depends on M^vi. From (fT7|) it is possible to lower M^i as well as the 
reheating temperature required by the thermal leptogenesis, though too small M^vi needs 
tiny Ai to avoid a strong washout. 

The cosmic-ray spectra of the anti-matter and photon also depends on the R-parity 
violation. In this letter, the flavor structure is determined by the neutrino Yukawa coupling 
and A*. By obtaining more cosmic-ray data in future such as FGST and AMS-02 and 
refining the astrophysical knowledge, we may distinguish the R-parity violating operators. 
Furthermore, since the R-parity violating operator can be embedded in a high-scale models 
such as SU(5) GUT and heterotic string models, combining with the future collider data we 
might be able to study physics in high-energy scale. 
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